A  i  '  1.  f  • 


MEMORANDUM 

RM-4179-NIH 

JULY  1964 


A  ■  ■ 

NUMERICAL  SOLUTIONS  OF 
FUNCTIONAL  EQUATIONS  BY  MEANS 
OF  LAPLACE  TRANSFORM-VII: 
THE  MELLIN  TRANSFORM 

Richard  Bellman 


Thi'  itixt-ol illation  wa*  Mi|>|M>rtf<l  in  |>arl  li\  I’liMir  llfallli  NT\iir  l{(•>■<•ar«'h  <«ranl 
NiimU-r  < i.  from  th<-  Dixi'ion  .if  (o-'ii-ral  Mt-iliial  S  ionci"-.  National 
Inotituti-s  of  (fcalth. 

dim:  awii.ahii.it^  notk  k 

Oualiticil  ri-ijiii'^ti'r-  max  olitain  ro|iii'>.  of  thi'  r)-|H>rt  from  tho  Dcfni'i'  IKh  nmi-ntation 
(a-ntfr  ( l)l)C). 


l>QO  ««IN  %f  •  lANV*  MONICA  •  CAII»0«NIA 


-itt- 


This  Memorandum,  with  Parts  I— VI,  comprises  a  study 
of  the  use  of  computing  machines  to  resolve  the  complex 
(.quations  which  arise  in  biology  and  medicine.  Part  VII 
discusses  the  question  of  the  numerical  inversion  of  the 
Mellin  transform. 
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SUMMARY 

In  this  series  of  papers,  we  have  so  far  considered 
various  analytic  and  computational  aspects  of  the 
numerical  inversion  of  the  Laplace  transform  and,  in 
addition,  some  applications  to  scientific  problems.  In 
this  paper  we  wish  to  consider  the  question  of  the 
numerical  inversion  of  the  Mellin  transform, 

s-1 

F(s)  ■  u(t)t  dt. 

0 

This  transform  also  plays  an  important  role  in  the  study 
of  functional  equations  of  the  type  that  arise  in 
mathematical  physics. 
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NUMERICAL  SOLUTIONS  OF  FUNCTIONAL  EQUATIONS  BY  MEANS  OF 
LAPLACE  TRANSFORM-^II:  THE  MELLIN  TRANSFORM 

1.  INTRODUCTION 

In  this  series  of  papers,  we  have  so  far  considered 
various  analytic  and  computational  aspects  of  the 
numerical  inversion  of  the  Laplace  transform  and,  in 
addition,  some  applications  to  scientific  problems.  In 
this  paper  we  wish  to  consider  the  question  of  the 
numerical  inversion  of  the  Mellin  transform, 

(l.I)  F(s)  -  u(t)t®  Mt. 

0 

This  transform  also  plays  an  important  role  in  the  study 
of  functional  equations  of  the  type  that  arise  in 
mathematical  physics. 

2.  A  PRELIMINARY  TRANSFORMATION 

Let  us  make  the  change  of  variable 


(2.1) 


t 

nrt' 


Then  (1.1)  assumes  the  form 


0 
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upon  simplifying  the  notation  by  setting 

(2.3)  g(r)  =  u(r/(l-r))/(l  -  r)^. 

We  suppose  that  u(t)  -*0  as  t  ->  or  rapidly  enough  so 
that  F(s)  is  defined  for  Re(s)  0,  and  thus  that 
g(l)  =  0. 

Applying  an  N— th  order  Gaussian  quadrature  as  in 
the  earlier  papers,  we  have 

N  r^.  s— 1 

(2.4)  F(s)  =r  w.g(r.) 

i«l  1-  1  i  -  r.. 

where  the  r^  are  the  roots  of  the  shifted  Legendre 
polynomial,  ^)  =  '/'j^(x). 

Setting  s  *  1,2,...,N,  we  obtain  N  linear 
algebraic  equations  for  the  N  quantities  w^g(r^), 
i  “  1,2,...,N,  Since  the  matrix  ([r^(l  —  r^)]-^)  is 
ill-conditioned,  this  is  not  necessarily  equivalent  to 
obtaining  meaningful  valuet  for  the  w^g(r^). 

3.  EXPLICIT  MATRIX  INVERSION 

Let  us  then  employ  a  device  used  in  [1]  to  obtain 
an  explicit  analytic  solution  of  the  corresponding  system 
obtained  in  connection  with  the  Laplace  transform.  From 

(2.4) ,  we  have 


(3.1)  Z  V,  -  Z  w^g(r.)['^' 

k-1  ^  ^  i-1  ^  ^  ^k-0  -' 

for  any  N  constants  a^,a2, » • • We  can  write  (3.1) 
in  the  form 


il  Vk  -  J, 


where 


k 

(3.3)  (p(x)  «  I  a,  ..x^, 

k-0 


a  polynomial  of  degree  N  —  1. 

Let  us  choose  this  polynomial  in  a  convenient 
fashion.  Suppose  that 


(3.4)  (p(  . 


Then,  if  y  ^  “  1,2,...,N,  are  the  corresponding 

coefficients  in  (3.3),  we  have 


N 

(3.5)  E  ^ki^k  “  Wig(ri)' 

k=i  J  J 

Thus,  (^+1  j)  required  inverse  matrix. 

To  obtain  the  required  polynomial  (pj(x),  consider 


the  function 


(3.6)  tTj^(x)  -  t 

where  ^j^(x)  is  the  function  introduced  in  Sec,  2.  We 
see  that 

(3.7)  i  - 

However,  ^j^(x)  is  a  rational  function,  not  a  polynomial. 
To  remedy  this,  consider  the  new  function 

(3.8)  Mj^(x)  -  (1  +  x)^tTj^(x)  -  (1  +  x)\^(^-]^  -), 


a  polynomial  of  degree  N,  Then 


(3.9)  <Pj(x) 

is  the  required 


I^(x) 


polynomial  of  degree  N  —  1. 


Writing 


(3.10) 


M 


N^l 

E 

k-0 
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we  obtain  the  desired  weighting  factors. 

4.  COMPUTATIONAL  ASPECTS 

A  computer  program  will  be  required,  as  in  the  case 
of  the  Laplace  transform  (see  [1]),  to  determine  the 


y  They  can,  in  this  way,  be  obtained  to  arbitrary 
precision.  We  shall  report  on  this,  together  with 
applications,  in  a  subsequent  paper. 
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